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Abstract 

A classical explanation of interference effects in the double slit experiment is proposed. We 
claim that for every single "particle" a thermal context can be defined, which reflects its em- 
bedding within boundary conditions as given by the totality of arrangements in an experimental 
apparatus. To account for this context, we introduce a "path excitation field", which derives from 
the thermodynamics of the zero-point vacuum and which represents all possible paths a "particle" 
can take via thermal path fiuctuations. The intensity distribution on a screen behind a double slit 
is calculated, as well as the corresponding trajectories and the probability density current. The 
trajectories are shown to obey a "no crossing" rule with respect to the central line, i.e., between 
the two slits and orthogonal to their connecting line. This agrees with the Bohmian interpretation, 
but appears here without the necessity of invoking the quantum potential. 
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1. INTRODUCTION 



Unless one assumed a purely pragmatic attitude and considered the formalism of quantum 
theory only as a set of rules to account for the empirical data, the challenge since the early 
days has been about developing a physical theory with a clear statement on what it is that 
is moving, and under which relevant boundary conditions and forces. 

Actually inherited from debates about the nature of light in still earlier centuries, this 
question was also in the twentieth century largely characterized by whether quantum sys- 
tems were, at the most fundamental level, particles or waves. A notable exception to these 
monistic approaches is the de Broglie-Bohm interpretation of the quantum mechanical for- 
malism with its model of particle plus "guiding" wave, i.e., a dualistic approach, albeit with 
an inherent asymmetry: the wave field acts on the particle, but not vice versa, a feature 
unknown to any other part of physics. 

Still, the modelling, or even the ontology, of all these approaches implies the existence of 
fundamental entities, be they particles or waves or something more complicated like, e.g., 
"wavicles", or "beables", etc.. These entities are assumed to characterize physical processes 
over a vast range of orders of magnitude, including the yet unknown domains of very small 
space and time scales down to the Planck length (where, eventually, some new physics may 
be required). However, these are unwarranted assumptions only. 

On the contrary, it may well be that there is practically an (almost?) infinite range of 
orders of magnitude, compared to the scales of today's highest resolutions. Facing such a 
scenario, or similar ones, it would be highly surprising, and even implausible, that at this 
point in history we should already know the "most fundamental level" of the laws of nature, 
i.e., in the form of the quantum physical laws. Note, for example, that if standard quantum 
field theory were valid up to the Planck scale, one would obtain the "famous" discrepancy 
of 120 orders of magnitude between theoretical and experimental values of the vacuum 
energy [1]. It would rather be much more plausible, then, that quantum phenomena, like all 
other physical phenomena we presently observe, refer to processes on some "intermediate 
levels". If the latter were the case, then "the quantum" would be an emergent phenomenon, 
i.e., based on some underlying sub-quantum physics, and limited by the boundary conditions 
of what one often terms "macroscopic physics" . 

In other words, like so many phenomena in the natural world, also the quantum phenom- 
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ena may result from the interplay between various levels of description, rather than being 
reducible to one "fundamental level" only. In this sense, we shall here try to model the 
quantum as a well-coordinated emergent system, where particle-like and wave-like phenom- 
ena are the result of both stochastic and regular dynamical processes. A prototype of such a 
system is well known from classical physics, viz., the "bouncing droplets" of Couder's group 
[2H6] , which in fact exhibit a whole series of phenomena reminiscent of quantum ones. Anal- 
ogously, our group has in recent years attempted to model a quantum as a nonequilibrium 
steady-state maintained by a permanent throughput of energy. Specifically, we consider 
a "particle" as a "bouncer" whose oscillations are phase-locked with those of the energy- 
momentum reservoir of the surrounding "vacuum", i.e., the zero-point field. (The possible 
existence of a corresponding, underlying "medium" is a priori independent of quantum the- 
ory. For similar views, compare, for example, the works of Hestenes [7j, Khrennikov [8], 
Nelson [9], and Nottale [10], as well as the approach of "stochastic electrodynamics" by 
de la Pena and Cetto, and in particular a recent paper on the "genesis of quantum nonlocal- 
ity" [llj.) In assuming that (part of) the zero-point fluctuations undergo regular oscillatory 
motion, where the latter is partly caused by and dynamically coupled to the oscillator's 
frequency u, we were able to derive fundamental elements of quantum theory from such a 
classical approach |1 2]. In ref. [13] we showed how Planck's relation between the energy E 
and the frequency u, E = hu, can be derived from our sub-quantum physics, with Planck's 
(reduced) constant h indicating a universal angular momentum, and we also showed that 
with this relation alone one can derive the exact Schrodinger equation from classical physics 
[lll[T5]. Moreover, also the stochastic element of the zero-point fluctuations enters decisively 
into our model, such that in effect we obtain an exact description of free quantum motion 
via a combination of the propagation of classical waves with diffusion due to stochastic sub- 
quantum mechanics. In particular, free quantum motion was shown to be exactly equivalent 
to anomalous ("ballistic") diffusion, which is a) crucially dependent on the initial conditions, 
but b) time-reversible despite the usual irreversible nature of the more common forms of 
diffusion processes. [16] 

To provide a suitable illustration, let us briefly discuss the following idealized scenario. 
Consider the straight path of a "particle" (e.g., an electron or a neutron) unaffected by 
any forces or apparatuses' boundary conditions. In other words, consider this path as a 
classical one. Then classical physics tells us that the Hamiltonian flow of such paths is 
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formally accompanied by Huygens waves originating from the particle, such that their wave 
fronts are always orthogonal to the particle's path. If, for example, these waves then enter 
a slit system, Huygens' principle tells us how the diffraction at the slit is to be described. 
Now assume that said Huygens waves are not just mathematical tools, but physical waves 
excited by a bouncer's oscillatory motion in its surrounding medium. Then, in this highly 
idealized classical scenario, Huygens' principle would describe the actual physical process 
of waves centrally-symmetrically originating from the "particle" position. As the bouncing 
oscillation is phase-locked with that of the waves, one can imagine a steering mechanism 
in analogy to that observed by Couder's group [SHH] such that the bouncer is actually a 
"walker" along its classical path. 

Now let us zoom in on such a system and consider the "particle" more closely, i.e., with 
statistical physics. If the bouncer were the only one in the world, with no obstacles or 
other boundaries affecting its movement, then the latter would - in its temporal mean - 
always remain on a straight path, with regular Huygens waves emanating from it. But in 
the real world there are other objects and boundary conditions, which will affect the above- 
mentioned medium surrounding the "particle", i.e., contribute to a physical "context", just 
like the walker's paths of Couder's group are affected by objects such as a slit system. This 
means for our approach that the waves emanating from the bouncer and those impinging on 
it from the surrounding world will interact: some, more regular, parts may create standing 
waves, e.g., due to a synchrony between incoming and outgoing waves, but many will just 
act as "disturbances" of the regular path, i.e., factually act as stochastic influences due to 
the irregular part of the zero-point fluctuations. In other words, Huygens' principle will 
then strictly apply only to the combination of a "walking" motion of some velocity v with 
centrally-symmetric diffusion waves' motion of velocity u orthogonal to it. It will thus be 
only approximate, disturbed by that part of the accompanying diffusive process which is to 
be described by the corresponding velocity fluctuations 6u. In fact, as we have shown for 
a Gaussian slit, the exact quantum mechanical result can be described as a combination of 
classical wave mechanics with the addition of a corresponding stochastic diffusion process 
[16] . In this paper, we intend to show that the same modelling procedure also perfectly 
applies to a double slit system, again providing the exact reproduction of the quantum 
mechanical results, albeit again relying on classical physics. 
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2. PATH EXCITATION FIELD AND DIFFRACTION AT A GAUSSIAN SLIT 



To begin with, we recall some of the basic results of our earlier work, including that on 
diffraction at a single Gaussian slit [H]. We claim for the system of a bouncer of frequency 
u embedded in a stochastic ( "zero-point" ) environment that its average total energy is given 
by the average "total" energy hu of the bouncer itself plus a kinetic energy term due to 
momentum changes =: mu which it receives from or gives off to the environment: 



Etot = ^u; + = const, (2.1) 
2m 

where the averaging (as denoted by the bars) is defined in n-dimensional configuration space 



as 



pl:= PplcTx. (2.2) 



P = P(x, t) refers to the probability density of some relevant distribution. For our model 
system, the latter will be introduced below as a solution of a generalized ("anomalous") 
diffusion equation, 

-P{x,t)=Dit)—P{x,t), (2.3) 
with D(t) being a time-dependent diffusion coefficient. As can be seen from ( 2.1[ ), the 



momentum changes can be either positive or negative, and actually will on the average 
balance each other, since they are a priori unbiased. The deeper reason for this balancing 
in our model is due to the fact that we consider the bouncer /walker to be a steady-state 
system in the sense that its "total" energy hw is maintained over times t ^ l/u hj the 
permanent throughput of kinetic energy p^/2m. In other words, to maintain the steady- 
state, during the intervals of the average order of t ^ there will both be an absorption 
of a momentum pu = mu and a release of the same amount, pu = —mu. 

In our earlier papers P^] - [TB] . or alternatively in [T7], for example, one derives from clas- 
sical physics the expression for u as 

VP 

u = -D^, (2.4) 

with the diffusion constant D. Noting that in reference [13], Planck's (reduced) constant h 
is given by a basic angular momentum, 

h = mcoT^, (2.5) 



and if r in (2.5) is identified with the usual expression for a diffusion length, r = y^2ZJ/u;, 



then one obtains the well-known diffusivity 

D = —. (2.6) 

2m ^ ^ 

So, apart from the ordinary particle current J(x, t) = P(x, i(:)v, we are dealing with two 
additional, yet opposing, currents J„, which are on average orthogonal to J P^HTB] . and 
which are the emergent outcome from the presence of numerous corresponding velocities 

^ VP 

",!, = "^^- < ' 

We denote with u+ and u_, respectively, the two opposing tendencies of the diffusion 
process. In the reference frame of a single free particle, and starting at t = at the center 
of the distribution P, the averages obey 

u_(x,t) = -u+(-x,t). (2.8) 

Now let us consider an experimental setup with a particle source. To describe the velocity 
distribution, we introduce a velocity field with average velocity v, and amplitudes -R(x, t). 



We refer to their intensities P = R as the solutions of a diffusion equation such as (2.3). 
These typically appear in the form of a Gaussian distribution P(x, t) of possible particle 
locations x, even if there is only one particle at a time emerging from the corresponding 
"Gaussian slit", in one dimension for simplicity, 

X (x — Xq)'^ 

P(a;,t) = — ^e-^^, (2.9) 



with the usual variance cr^ = (Ax)^ = (x — Xq)^, and where we choose xq = further on. 
Regarding u, even in this scenario of one-particle-at-a-time, we deal with an ensemble of 
velocity vectors Ua{t) representing hypothetical motions on the sub-quantum level in a small 
volume around x, whose mean value will be given by 

^ 7V(x,t) 

"(x,t) = ^^ E ""W- (2.10) 

Here, the (typically very large) number refers to the number of possible path directions of 
the bouncer due to the existence of the wave-like excitations of the zero-point field. This may 
be reminiscent of Feynman's picture of photons virtually probing every possible path in an 



experimental setup, but in our case it is a configuration of real wave-like excitations, which 
in a resulting Brownian-type motion guide the bouncer along its path of average velocity 
u. Again, we note that here we discuss not only a passive guidance of the "particle" by the 
surrounding wave configurations, but point out also the very active role of the "particle" in 
(partly) creating said wave configurations due to the effects of its bouncing. To account for 
( 2.7[ ), we split up u(x, t) according to 

N N 



U X, t = 



^[u+ + u_], (2.11) 



thus reflecting the isotropy of the diffusion process. Still, the uncontrollable and possibly 
unknowable velocity field u representing the Brownian motion of the bouncer may not 



be operational, but when we take the average according to the rule (2.2), we obtain a 
"smoothed-out" average velocity field 



u(x, t) = y Pu(x, t) d^x, (2.12) 

which is all that we need for our further considerations. Similarly, based on the fact that we 
have an initial Gaussian distribution of velocity vectors v(x, t), we define an average velocity 
field V of the wave propagation as 

v(x,t) = J Pv(x,t)d"a;, (2.13) 

and we here just note what we have discussed extensively in previous papers [TIHTUI [TH] . 
i.e., that there is an average orthogonality between the two velocity fields, u and v, 

= y Pv ■ u d"x = 0. (2.14) 

This holds for a variety of reasons, one of them simply being that u and v are linearly 
uncorrelated. 

In effect, then, the combined presence of both velocity fields u and v can be denoted as a 
path excitation field: via diffusion, the bouncer in its interaction with already existing wave- 
like excitations of the environment creates a "heated-up" thermal "landscape", which can 
also be pictured by interacting wave configurations all along between source and detector 
of an experimental setup. Recall that our prototype of a "walking bouncer", i.e., from the 
experiments of Couder's group, is always driven by its interactions with a superposition of 



waves emitted at the points it visited in the past. Couder et ah denote this superposition 
of in-phase waves the "path memory" of the bouncer IE]. This imphes, however, that the 
bouncers at the points visited in "the present" necessarily create new wave configurations 
which will form the basis of a "path memory" in the future. In other words, the wave 
configurations of the past determine the bouncer's path in the present, whereas its bounces 
in the present co-determine the wave configurations at any of the possible locations it will 
visit in the future. Therefore, we call the latter configurations the path excitation field, which 
may also be described as "heated-up" thermal field. As in the coupling of an oscillator with 
classical diffusion, diffusion wave fields arise with instantaneous field propagation [151 US] , 
one has elements of the whole setup which may be nonlocally oscillating ("breathing") in 



phase. This means that the Gaussian solving (2.3) does represent a nonlocal path excitation 



field in that it is a physically existing and effective entity responsible for where the bouncing 
"particle" can possibly go. As we have shown [16], one can in this classical framework, along 
with the time-dependence of D(t), effectively and easily describe the (sub-) quantum physics 
of diffraction at a Gaussian slit, which we now briefly recapitulate. 



At first we note that Eq. (2.1) is an average energy conservation law only. This means 
that apart from the momentum changes = ±mu discussed so far, also variations in pu 
will have to be taken into account, and thus also variations in the "particle energy" hu. If 
for the latter one just considers its kinetic energy term, mv'^/2, then said variations will lead 
to exchanges of velocity/momentum terms providing the net balance 

m5v = mSu. (2.15) 

Thus, the average kinetic energy variation will be determined, once an initial condition is 
defined: 



5Ekin(t) = 5^kin(0) = const. (2.16) 

In other words, one starts with the initial average squared momentum change, which is 
determined by the characteristics of the particle source (to be shown below), 

=: m\l, (2.17) 



2 2 2 

pf^ = m 



t=o 



relative to which all further momentum balances have to be carried out. Then, (2.16) reads 
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Using the expression (2.7) for u, one obtains with Eqs. (2.9) and (2.6), and with (VlnP)^ 



- V2 In P, that 



{6u) 



a" 



+ {6uy 



(2.19) 



where as usual cr = a{t) = y for xo(t = 0) = 0, and ctq = ^"(t = 0). One can view the 
Gaussian distribution P of kinetic energy also as a sort of "heat accumulation" , which has 
its maximum at the center. Considering now the application of momentum fluctuations (up 
to second order) to a particle with initial (t = 0) distance x(0) from said center, with the 
fluctuation term for t > defined as pu ± 5pu = ±m(u ± 5u), one will find said particle at 
time t at the location 

x(t) = x(0) ± (u ± (5u) t. (2.20) 



Squaring Eq. (2.20) and forming the r.m.s. then provides with Eq. (2.18) and average or- 



thogonality (2.14), and thus also u ■ x(0) = [TU], 

^2(t) = ^(0) + + (5^1 = ^2(0) + ult"^. 



(2.21) 



(Note that we could have included still higher order terms of the fluctuations, but due to 



the conservation condition (2.16), Eq. (2.21) results also for any number of higher orders.) 



Comparing with Eq. (2.19) also provides the time evolution of the wave packet's variance 
as 



a 



cr, 







(2.22) 



As the spreading ratio aja^ for the wave packet is independent of x, and thus not only valid 



for the particular point x{t) = cr(t), one can generalize the square root of Eq. (2.22) for all 
x, thus providing 

x{t) = x(0) — , where 



a 



1 + 



— ■ (2-23) 

CTq (To y (Tg 

Thus, one obtains for "smoothed-out" trajectories (i.e., averaged over a very large number 
of Brownian motions) a sum over a deterministic and a fluctuations term, respectively, i.e., 



with Eq. (2.6), 



a;tot(t) =vt + x{t) =vt + x{0)— =vt + x(0)Wl + 



To 



(2.24) 



Moreover, one thus derives with Eq. (2.19) the average velocity field of a Gaussian wave 
packet as 

Vto,{t) = v{t) + = v{t) + [xtotit) - vt] %. (2.25) 



dt 



Note that Eqs. (2.22)-(2.25) are derived solely from statistical physics. Still, they are 
in full accordance with quantum theory, and in particular with Bohmian trajectories 



Note also that one can rewrite Eq. (2.21) such that it appears like a linear-in-time formula 
for Brownian motion, 

= x^{0) + D{t) t, (2.26) 

where a time dependent diffusivity 



D{t) = u^t 



t 



(2.27) 



characterizes Eq. (2.26) as ballistic diffusion. The appearance of a time-dependent D{t) is 
easy to understand. The diffusivity is changed over time, because the "particle's" thermal 
environment changes: With the "heat" initially concentrated within the narrow spatial con- 
straints determined by ctq of the source ("Gaussian slit"), D{t) must become larger with time 
because of the gradually lower heat concentration due to dissipation into the unconstrained 
environment. (Note that a similar scenario was suggested by Garbaczewski ^17j; others are 
presently intensively discussed in the context of the so-called "superstatistics" [2T]. More- 
over, ballistic diffusion has also been measured recently in experiments realizing quantum 
walks, see, e.g., |21| and references therein.) This makes it possible to simulate the disper- 
sion of a Gaussian wave packet on a computer by simply employing coupled map lattices 



for classical diffusion, with the diffusivity given by Eq. (2.27). (For detailed discussions, see 
refs. [IS] and [ig.) 



3. AN EXPLANATION OF INTERFERENCE EFFECTS IN THE DOUBLE SLIT 
EXPERIMENT 

With the essentials of Gaussian dispersion at our disposal, it is straightforward to now 
also describe and explain quantum interference with our approach. We choose a textbook 
scenario in the form of the calculation of the intensity distribution and the particle trajec- 
tories in an electron interferometer. As we are also interested in the trajectories, we refer 
to, and compare our results with, the well-known work by Philippidis et al. [23], albeit in 
the form as presented by Holland [20] . 

We choose similar initial situations as in [20], i.e., electrons (represented by plane waves 
in the forward ?/-direction) from a source passing through "soft-edged" slits 1 and 2 in a 
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barrier (located along the x-axis) and recorded at a screen. In our model, we therefore note 
two Gaussians representing the totality of the effectively "heated-up" path excitation field, 
one for slit 1 and one for slit 2, whose centers have the distances +X and —X from the 
plane spanned by the source and the center of the barrier along the y-axis, respectively. We 
recall from classical wave mechanics that the plane-wave superposition of two amplitudes 
Ri = y/Pi, i = 1 or 2, with resulting amplitude R and wave number k provides 



Rk = i?iki + i?2k2. 



(3.1) 



such that for |ki| = |k2| = |k|, one obtains (with hats further on denoting unit vectors) the 
total intensity 

« 2 



tot •- 



R' 



Riki + i?2k2 



Rl + Rl + 2RiR2Cosip = P^ + P2 + 2y^PiP2Cosip. (3.2) 



The x-components of the centroids' motions from the two alternative slits 1 and 2, respec- 
tively, are given by the "particle" velocity components 



m 



(3.3) 



respectively, such that the relative group velocity of the Gaussians spreading into each other 
is given by Av^ = 2vx. However, in order to calculate the phase difference (f descriptive 



of the interference term of the intensity distribution (3.2), one must take into account also 



the wave packet dispersion as described in the previous Chapter. Thus, one obtains with 



Eq. (2.25) the total relative velocity of the two Gaussians as 

t 



ot,x 



-{X + v^t) 



Aim? (T^cTq 



(3.4) 



Therefore, the total phase difference between the two possible paths (i.e., through either 
slit) becomes 

where the last term also reads as 



^(mAt;tot,x x) = 2mv^^ - {X + ^^•^)^^^^' 



(3.5) 



-(X + Vxt)x 



1 ult 
D 



-(X + Vxt)x 



1 a 
Da' 



The Gaussians Pi and P2 for the corresponding slits are given as 

1 



V27ia^ 
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[x-{X+v^t)]^/2a^ 



(3.6) 



(3.7) 



and 



[x+{X+v^t)]^/2a^ 



V2 



(3.8) 



With equal amplitudes R = Ri = y/Pi, for i = 1,2, of the Gaussians, and with normalization 
constant A^, we thus obtain the usual interference pattern in the form of the intensity 
distribution: 

1 



Ptot{x,t) 



R^N^ 



V2n 



+ 2cos¥?|, (3.9) 



with Lf given by Eq. (3.5). The term inside the curly brackets of Eq. (3.9) describes the 



interference fringes whose "dark" nodes are at the locations well-known from textbooks 

1 



X 



(n + -)TT/k^, 



(3.10) 



if, according to Eq. (3.6), either a/a is negligible (i.e., a{t) ~ const.), or under the specific 
time- conditions X = —Vxt, as the wave packets must approach each other [vx < 0) for t > 0. 
For a discussion of further aspects including, for example, those on the Fraunhofer limit, we 
refer the reader to the details provided in Holland's book [20j, which equally result from our 
Eq. (Q. 

Fig. [T] depicts the interference of two beams emerging from Gaussian slits established with 
the aid of a simulation as in ref. [16], where we simulate diffusion with a time-dependent 
diffusivity D(t). To account for interference, we simply follow the classical rule for the 



intensities (3.2), with (p from Eq. (3.5). The trajectories are the flux lines obtained by 



choosing a set of equidistant initial points at y = 0. Two adjacent flux lines thereby define 
regions of constant difference AP. 

Just as with the dispersion of a single Gaussian, we want to stress that also the quantum 



interference pattern (3.9) has been derived here without the use of quantum mechanics, but 



solely on the basis of classical physics as it was outlined in the footnote of the introduction. 
Exploiting our concept of the "path excitation field" , we have for this derivation implicitly 
used the assumption of velocity fields u(x, t) and v(x, t), respectively, which have entered 



the expression of a Gaussian's average velocity field, Eq. (2.25), and which in turn were 



shown to essentially contribute to the interference pattern (3.9). In what follows, we now 



want to make the use of these velocity fields more explicit, i.e., we shall now concentrate on 
understanding the emerging particle trajectories during quantum interference on the basis 
of our classical velocity fields. 
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FIG. 1. Classical computer simulation of the interference pattern: intensity distribution with 
increasing intensity from white through yellow and orange, with trajectories (red) for two Gaussian 
slits, and with small dispersion (evolution from bottom to top; Vx,i = —fx, 2)- The trajectories follow 
a "no crossing" rule: particles from the left slit stay on the left side and vice versa for the right 
slit. This feature is explained here by a sub-quantum build-up of kinetic (heat) energy acting as 
an emergent repellor along the symmetry line. 

4. GEOMETRIC MEANING OF THE PATH EXCITATION FIELD AND ITS 
TRANSLATION INTO QUANTUM LANGUAGE 

Let us consider a single, classical "particle" ("bouncer") following the propagation of a 
set of waves of equal amplitude R, each representing one of i possible alternatives according 
to our principle of path excitation. We now concentrate on the specific role of the velocity 



fields u, which were present in Chapter [3] only implicitly in the expressions (3.2) and (3.5) 
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To describe the required details, each path i be occupied by a Gaussian wave packet with 
a "forward" momentum pj = Kki = m\i. Moreover, due to the stochastic process of path 
excitation, the latter has to be represented also by a large number of consecutive Brownian 



shifts, Y^u,a = rnua- Recalling (2.11), one has 

N 



u(x,t) 



1 

2N 



N 



U I + u 



(4.1) 



_a=l a=l 

This provides for the case of interference at a double slit the total velocity field (with indices 
i = 1 or 2 referring to the two slits) 



Vtot — Vtot,i + Vtot,2 Vi 



Ul 



Ul 



V2 



U2+ , U2 



2 2 2 2 

With two Gaussian distributions Pi = R\ and P2 = R2 as given above, one has 

-RtotVtot = -RlVtot.l + -R2Vtot,2 5 



(4.2) 



(4.3) 



which approximates the addition of amplitudes at the crossing point of plane waves. To ob- 
tain an expression for Ptot = -Rtot ^^e relevant velocity components in the interference 



region, we derive from Eq. (4.3) 

-Rtot = -Ri Vi + 



Ul+ , Ul 



ry I - U2+ U2- 



Vtot 



(4.4) 



2 2 

Taking into account the conservation of the "particle momentum" in both channels, we have 
I Vtot I = |vtot,i| = |vtot,2|, leading to 



P 



tot 



i?l V 



Ul- 



Ui 



+ i?2 V2 



U2H 



U2- 



Vtot- 



, , (4.5) 
2 2 J V 2 2 ' ^ ^ 

To help with the bookkeeping. Fig. |2] displays all the relevant vectors and some of their 



included angles. With Eq. (2.8), and setting ui+ ui and ui_ — t- — ui, one obtains the 
total average current 

Jtot = -R^ot^tot = -RiVi + -Riv2+ 

+R1R2' 



Vi + V2) COS(Vi, V2 



+ I Y + V2 ) COS(Ui, V2) 



Vl + 



Ul 



U2 



Ul _ U2 

2 2 



COSl^Ui, U2J 



Ul U2 

2 2 



cos(vi 


U2) 






U2 

2 


cos(ui 


V2) 


- ( 


^Ui 

v2 




cos(ui 


U2) 




Ul 

' 2 


U2 

2 



COS(^Vi,U2j 
C0S(U1,U2) 
COs(Ui,U2) f, 

(4.6) 
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FIG. 2. Two-dimensional schematic of the path excitation field behind a double slit. Uj and Vi, 
z = 1 or 2, denote the average directions of the velocity fields involved, with Vi referring to single free 
"particle" velocities (i.e., either through slit 1 or slit 2), and the Ui referring to additional diffusion 
velocities. The red arrows indicate how the mixing of velocity fields Vi and Uj, with i ^ j, produces 



the terms proportional to sin 99 in Eq. (4.7). Upon superposition of the two channels, according to 



Eq. (4.7), one obtains the final particle trajectories whose emergent behavior is characterized by a 



"no crossing" rule (see Fig. [T]). 



and thus finally 



PiVi + P2V2 + V P1P2 (vi + V2) cos V3 + V (ui - U2) sin (p. 



(4.7) 



Note that Eq. (4.7), upon the identification of Uj 



^ from Eq. (2.7) and with Pj = R 



turns out to be in perfect agreement with a comparable "Bohmian" derivation [201 El] • The 



formula for the particle trajectories, then, simply results from Eq. (4.6), i.e. 

1 ^ 



P 



•tot , 



(4.^ 



tot 



with Ptot given by Eq.(3.2) and J^t given by Eq. (4.7). 
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In Fig. [T| one can observe a basic characteristic of the (averaged) particle trajectories, 
which, just because of the averaging, are identical with the Bohmian trajectories. To fully 
appreciate this surprising characteristic, we remind the reader of the severe criticism of 
Bohmian trajectories as put forward by Scully and others (see [25], and references therein.) 
The critics claimed that Bohmian trajectories would have to be described as "surreal" ones 
because of their apparent violation of momentum conservation. In fact, due to the "no 
crossing" rule for Bohmian trajectories in Young's double slit experiment, for example, the 
particles coming from, say, the right slit (and expected at the left part of the screen if 
momentum conservation should hold on the corresponding macro-level) actually arrive at 
the right part of the screen (and vice versa for the other slit). In Bohmian theory, this "no 
crossing" rule is due to the action of the non-classical quantum potential, such that, once 
the existence of a quantum potential is accepted, no contradiction arises and the trajectories 
may be considered "real" instead of "surreal" . 

Here we can note that in our sub-quantum approach an explanation of the "no crossing" 
rule is even more straightforward and actually a consequence of a detailed microscopic mo- 
mentum conservation. As can be seen in Fig.[T| the (Bohmian) trajectories are repelled from 
the central symmetry line. However, in our case this is only implicitly due to a "quantum 
potential" , but actually due to the identification of the latter with a kinetic (rather than a 
potential) energy: As has already been stressed in [TJ], it is the "heat of the compressed 
vacuum" that accumulates along said symmetry line (i.e., as reservoir of "outward" oriented 
kinetic energy) and therefore repels the trajectories. Fig. [T] is in full concordance with the 
Bohmian interpretation (see, for example, for comparison). However, as mentioned, in 
our case also a "micro-causal" explanation is provided, which brings the whole process into 
perfect agreement with momentum conservation on a more "microscopic" level. 

This can be shown even in greater detail. Whereas in the example of Fig. [T] the small 
amount of dispersion is practically negligible, we now discuss in more detail an interference 
scenario with significant dispersion of the two Gaussians, i.e., where initially the two Gaus- 
sians spread independently from each other due to the action of the diffusive velocities ui 
and U2, respectively. In Fig. [sj trajectories according to Eqs. (4.7) and (4.8) for the two 



Gaussian slits are shown. The interference hyperbolas for the maxima characterize the re- 
gions where the phase difference (p = 2n7r, and those with the minima lie at (p = {2n + 1)it, 
n = 0,±1,±2, ... Note in particular the "kinks" of trajectories moving from the center- 
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FIG. 3. Classical computer simulation of the interference pattern: intensity distribution with 
increasing intensity from white through yellow and orange, with trajectories (red) for two Gaussian 
slits, and with large dispersion (evolution from bottom to top; Vx,i = Vx,2 = 0). The interference 
hyperbolas for the maxima characterize the regions where the phase difference (p = 2mr, and 
those with the minima lie at (p = {2n + l)7r, n = 0, ±1, ±2, . . . Note in particular the "kinks" of 
trajectories moving from the center-oriented side of one relative maximum to cross over to join 
more central (relative) maxima. In our classical explanation of double slit interference, a detailed 
"micro- causal" account of the corresponding kinematics can be given. 

oriented side of one relative maximum to cross over to join more central (relative) maxima. 
The trajectories are in full accordance with those obtained from the Bohmian approach, as 
can be seen by comparison with references [20], [26] and [27], for example. In our classical ex- 
planation of double slit interference, a detailed "micro-causal" account of the corresponding 
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kinematics can be given: Firstly, we note that the last term in Eq. (4.7), which is responsi- 
ble for the genuinely "quantum" behavior, is characterized by the vector ui — U2 which in 
the interference region is always oriented into the "forward" direction away from the slits 
(Fig. [2]). This means that said last term is modulated by simp, with the result that the term 
alternates between "forward" directions where simp > and "backward" directions where 
sin p> < 0. 

Thus, in the "backward" cases, the trajectories coming from the relative maxima (bright 
fringes) loose velocity/momentum in the forward direction and cross over into the area of 
the relative minimum (dark fringes). Alternatively, in the "forward" cases, the trajectories 
coming from the relative minima (dark fringes) gain velocity /momentum in the forward 
direction and thus align with the other trajectories of the bright fringes. In other words, for 
the areas where simp < 0, part of the current (along a relative maximum) is being removed 
("depletion"), whereas for simp > 0, parts of currents flow together to produce a newly 
formed bright fringe ("accumulation"). This is in accordance with our earlier description of 
quantum interference, where the effects of diffusion wave fields were explicitly described by 
alternating zones of heat accumulation or depletion, respectively [15J. Towards the central 
symmetry line, then, one observes heat accumulation from both sides, and due to big mo- 
mentum kicks from the central accumulation of heat energy, the forward particle velocities' 
directions align parallel to the symmetry axis. With the crossing-over of particle trajecto- 



ries being governed by the last, diffusion-related, term on the right hand side of Eq. (4.7), 
one finds that for ip = the resulting diffusive current is zero and thus, as total result of 
the overall kinematics, no crossing is possible. Further, we note that our results are also 
in agreement with the recently published experimental results by Kocsis et al. [28j. Here 
we just comment that, as opposed to the Bohmian interpretation, we give a micro-causal 
explanation of these results without invoking the quantum potential. 

Finally, to illustrate the straightforward applicability of our model to more general situa- 
tions, i.e., as compared to the simple symmetrical scenarios of Figs. [T] and [3| we now employ 
our simulation schema to cases where neither the Gaussians are identical nor their weights. 
We thus study asymmetric coherent superpositions as discussed in reference |2l], and in our 
Figs. [4] to [6] we show results in accordance with the figures 4-6 of ref. [21]. The analysis of 
ref. [21] holds identically in our approach, so that we here restrict ourselves to pointing out 
that our figures display the following cases of varied properties for the beams emerging from 
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the two slits: 



• different modulus of the initial velocity/momentum, 

• different initial spreading, 

• different weights for the probability densities. 




FIG. 4. Same as Fig. [T| but with different initial average velocities: Vx,2 = —^Vx^i- Note again 
the "no crossing" behavior, with the two trajectory bundles repehing each other due to the kinetic 
(heat) energy along the slanted central line. 



5. THE SCHRODINGER EQUATION REVISITED 

Upon employment of the Madelung transformation for each path j [j = 1 or 2), 

^j. = i^e*^^/^ (5.1) 



and thus Pj = i?| = l^^^p = with the definitions ([2^ and Vj := VSj/m, ip 



{Si — S2)/h, and recalling the usual trigonometric identities such as cosy^ = 1 (e**^ + e "''), 
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FIG. 5. Same as Fig. |4j but with different initial spreading: cji = 3a2- Although the two par- 
tial beams altogether reflect off each other, one can clearly observe the effect of "microscopic" 
momentum conservation: the path excitation field of the right beam is propagated over to the 
micro-kinematics of the left beam, and vice versa. 




FIG. 6. Same as Fig. [5| but with different probability densities: Pi = 2P2. Note that the emerging 
beam behavior compares more with inelastic scattering than with the elastic-type scattering of 
Fig. [5} as part of the left beam merges with the right one. 
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one can rewrite the total average current (4.7) immediately as 



'tot 



-PtotVt 



ot 
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.,V(*i + ^2)A , 1 
-m^^ 1 + 



ih 
2m 



[\[/*V\l/ - vl/V\l/* 



m 



-Re{^*(-iW)^}, 



m 



ih 



V(^i + ^2) 



(5.2) 



where Pt 



tot 



1^1 + ^2! 



|\l'p. In the Appendix we show an alternative derivation of 



Eq. (5.2) from our classical approach. The last two expressions of (5.2) are the exact well- 



known formulations of the quantum mechanical probability current, here obtained without 



any quantum mechanics, but just by a re-formulation of (4.7). In fact, it is a simple exercise 



to insert the wave functions (5.1) into (5.2) to reobtain (4.7). 



Moreover, one obtains with (5.2) for the total average current the conservation law 
dPtot 



dt 



-V - J 



tot 



ih 

2m 



(5.3) 



As it also holds that 



dt 



dt 



(\]/*vI>) = vl>*vi> -|- vl/vj/* 



a comparison of Eq. (5.3) and (5.4) provides 



ih 



v]/*v]/ -|- v]/v]>* 



h^ 

2m 



(5.4) 



(5.5) 



which is the difference between two equations, i.e 

^2 



and 



2m 



h^ 

V^^* + V"^* 1 ^ 

2m 



-ih^!*^! = 

and thus of the Schrodinger equation (with external potential V) 

h^ 

ihm = v^^ + 

2m 



(5.6) 



(5.7) 



(5.8) 



and its charge conjugate. 

The validity of the Schrodinger equation is here seen to be based on the assumption 
constitutive for our model, i.e., of the quantum being an off-equilibrium steady-state system, 
with its oscillations cu being maintained by a steady and balanced throughput of energy. 
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This raises, of course, the possibihty to consider the breakdown of that balance, which is 
conceivable at time scales t c^i^l/u oi shorter. For, we have assumed throughout that said 
balance is maintained between diffusive currents, which are responsible for the maintenance 
of the oscillation of frequency u. So, an eventual breakdown of the Schrodinger equation and 
the usual laws of quantum mechanics can be assumed when the symmetry between diffusive 
currents is broken, i.e., on the mentioned timescales. 

6. CONCLUSIONS AND OUTLOOK 

Considering "particles" as oscillators ( "bouncers" ) coupling to regular oscillations of the 
"vacuum's" zero-point field, which they also generate, we have shown how a quantum can 
be understood as an emergent system. In particular, the dynamics between the oscillator 
and the "bath" of its thermal environment can be made responsible not only for Gaussian 
diffraction at a single slit [16], but also for the well-known interference effects at a double 
slit. We have also shown how the model entails the existence of a path excitation field, i.e., a 
field spanned by the average velocity fields v(x, t) and u(x, t), respectively, where the latter 
refer to diffusion processes reflecting also the stochastic parts of the zero-point field. We 
have derived in this paper, on the basis of classical physics, i.e., without the use of quantum 
mechanical terms, the exact intensity distribution at a screen behind a double slit, as well as 
the details of the more complicated particle current, or of the Bohmian particle trajectories, 
respectively. 

In a computer simulation, moreover, we have pictured quantum interference employing 
well-known techniques to model classical diffusion processes. With this new tool at hand, 
and having proven the identity with the results of the usual quantum mechanical formalism, 
we are encouraged to tackle further foundational issues, such as relativistic formulations, 
and, in particular, questions concerning nonlocality and entanglement. As to the latter, one 
thing can be noted right away. Considering double interferometry, for example, with the 
emission of two anti-correlated particles with opposite directions, the nature of the path 
excitation field immediately implies two features: a) the "wholeness" of the processes going 
on in the total domain of the experimental setup, and b) nonlocal correlations between 
distant parts of the experiment. 

For, it is obvious that, say, the very existence of the initial Gaussian representing the 
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source (and at some time "decaying" into two Gaussians propagating in opposite directions), 
implies one path excitation field in both directions. This means that whereas the average 
velocities for the respective particles, v and v' are anti-parallel and refer to separate regions 
in space, it is one and the same diffusive velocity field u = u' that affects both particles 
throughout the experimental setup! We shall show in a forthcoming paper how this feature 
of "wholeness" implies the existence of nonlocal correlations. Due to the nonlocal nature of 
the involved diffusion wave fields, and based on our proposed model, it should be possible 
to prove a corresponding identity with entangled states in quantum mechanics. 

The phenomenon of entanglement is thus possibly rooted in the existence of the path 
excitation field, a version of which we already encountered in the present paper, i.e., with 
the usual double slit interferometry. So, one can expect that the concepts developed here 
will also help to deepen our understanding of entanglement. As far as explanations of the 
physics of the double slit experiment are concerned, however, we think that we can safely 
say that it is possible to understand quantum mechanics. 



Appendix A 



In this appendix, we give an alternative derivation of the quantum mechanical current 



(5.2) from classical physics. We have in previous papers employed the Madelung transfor- 



mation 



(Al) 



and thus P = E? = |\E'p = \E'*\E', to obtain a "vocabulary" for the translation of our notation 
into that of quantum mechanics. 
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which provides also that [5^ 
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(A2) 



(A3) 



(Note that generally the division by \E' is regular, since wherever u is defined, a particle path 
exists, and therefore 7^ 0.) Thus, one can also write down the momentum vector p^.^^ as 
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a complex one providing 
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tot 



(In quantum mechanical terms, Eq. (A4) is the corollary of the fact that p^o^ is the "eigen- 
vector" of the "momentum operator" p = —iKV.) Moreover, as was shown in [29], one can 
use the classical wave mechanics formula (3.1), -Rktot = -Riktot,i 



-R2ktot,2, and Eq. (A2) for 



the case of two alternative paths 1 and 2, to write k'^t 



/Pktofktotd' 
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corresponding to ^tot = + ^E^s 



If, with these ingredients, one now defines the average velocity Vtot 
immediately has 

Vtot = 
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from which one reobtains the total quantum mechanical current (5.2). 
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